Early-universe constraints on Dark Energy 
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In the past years 'quintessence' models have been considered which can produce the accel- 
erated expansion in the universe suggested by recent astronomical observations. One of the 
key differences between quintessence and a cosmological constant is that the energy density in 
quintessence, £7^, could be a significant fraction of the overall energy even in the early uni- 
verse, while the cosmological constant will be dynamically relevant only at late times. We 
use standard Big Bang Nucleosynthesis and the observed abundances of primordial nuclides 
to put constraints on fi^ at temperatures near T IMeV . We point out that current ex- 
perimental data does not support the presence of such a field, providing the strong constraint 
r20(MeV) < 0.045 at 2cr C.L. and strengthening previous results. We also consider the effect 
a scaling field has on CMB anisotropies using the recent data from Boomerang and DASI, 
providing the CMB constraint fl^ < 0.39 at 2a during the radiation dominated epoch. 



Introduction. Recent astronomical observations [|l| 
suggest that the energy density of the universe is 
dominated by a dark energy component with nega- 
tive pressure which causes the expansion rate of the 
universe to accelerate. One of the main goals for cos- 
mology, and for fundamental physics, is ascertaining 
the nature of the dark energy |^ . 

In the past years scaling fields have been con- 
sidered which can produce an accelerated expansion 
in the present epoch. The scaling field is known 
as "quintessence" and a vast category of "tracker" 
quintessence models have been created (see for exam- 
ple and references therein), in which the field 
approaches an attractor solution at early times, with 
its energy density scaling as a fraction of the domi- 
nant component. The desired late time accelerated 
expansion behaviour is then set up independently of 
initial conditions, with the quintessential field domi- 
nating the energy content. 

Let us remind the reader of the two key differ- 
ences between the general quintessential model and 
a cosmological constant: firstly, for quintessence, the 
equation-of-state parameter = p/ p varies in time, 
usually approaching a present value < —1/3, 
whilst for the cosmological constant remains fixed at 
w\ = — 1. Secondly, during the attractor regime the 
energy density in quintessence Q.^ is, in general, a sig- 
nificant fraction of the dominant component whilst JIa 
is only comparable to it at late times. 

Future supernovae luminosity distance data, as 
might be obtained by the proposed SNAP satel- 
lite, will probably have the potential to discriminate 
between different dark energy theories These 
datasets will only be able to probe the late time be- 
haviour of the dark energy component at red-shift 
z < 2, however. Furthermore, since the luminosity 
distance depends on w through a multiple integral re- 
lation, it will be difficult to infer a precise measure- 
ment of w{z) from these datasets alone 

In this Letter we take a different approach to the 
problem, focusing our attention on the early time be- 
haviour of the quintessence field, when the tracking 



regime is maintained in a wide class of models, and 
ri^ is a significant (> 0.01, say) fraction of the overall 
density. 

In particular, we will use standard big bang nucle- 
osynthesis and the observed abundances of primordial 
nuclides to put constraints on the amplitude of fJ^ 
at temperatures near T ~ IMeV. The inclusion of a 
scaling field increases the expansion rate of the uni- 
verse, and changes the ratio of neutrons to protons 
at freeze-out and hence the predicted abundances of 
light elements. 

The presence of this field in the radiation domi- 
nated regime also has important effects on the shape 
of the spectrum of the cosmic microwave background 
anisotropies. We use the recent anisotropy power 
spectrum data obtained by the Boomerang [Q and 
DASI 1^ experiments to obtain further, independent 
constraints on U,^ during the radiation dominated 
epoch. 

There are a wide variety of quintessential models; 
we limit our analysis to the most general ones, with 
attractor solutions established well before nucleosyn- 
thesis. 

More specifically, we study a tracker model based 
on the exponential potential V — Vbe^'^'^ H. If the 
dominant component scales as p„ — /5o(^"i then 
the scaling field eventually approaches an attractor 
solution, and its fractional energy density is given by 
^4> = However, the pure exponential potential, 
since it simply mimics the scaling of the dominant 
matter in the attractor regime, cannot produce an ac- 
celerated expanding universe in the matter dominated 
epoch. 

Therefore, we focus our attention on a recently pro- 
posed model by Albrecht and Skordis (referred to as 
the AS model from herein) , motivated by physics 
in the low-energy limit of M-theory, which includes 
a factor in front of the exponential, so that it takes 



the form V ^ ((0o - + A) , 



-A0 



The prefactor 



introduces a small minimum in the potential. When 
the potential gets trapped in this minimum its kinetic 
energy disappears, triggering a period of accelerated 
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FIG. 1. Top panel: Time behaviour of the fractional 
energy density for the Albrecht and Skordis model to- 
gether with the constraints presented in the paper. The 
parameters of the models are (assuming h = 0.65 and 
= 0.65) A = 3, (^0 = 87.089, A = 0.01 and A = 8, 
00 = 25.797, A = 0.01. Bottom panel: Time behaviour 
for the overall equation of state parameter wtot for the 
two models. Luminosity distance data will not be useful 
in differentiating the two models. 



expansion, which never ends if < 1 . 

In Fig. H we introduce and summarise the main re- 
sults of the paper. In the figure, the BBN constraints 
obtained in section 2, and the CMB constrains ob- 
tained in section 3 are shown together with two dif- 
ferent versions of the AS model which both satisfy the 
condition fl^ = 0.65 today. 

Constraints from BBN. In the last few years impor- 
tant experimental progress has been made in the mea- 
surement of light element primordial abundances. For 
the '^He mass fraction, ^jjg, two marginally compati- 
ble measurements have been obtained from regression 
against zero metallicity in blue compact galaxies. A 
low value Yfjg = 0.234 ± 0.003 [||| and a high one 
Yjjg = 0.244 ± 0.002 give realistic bounds. We 
use the high value in our analysis; if one instead con- 
sidered the low value, the bounds obtained would be 
even stronger. 

Observations in different quasar absorption line sys- 
tems give a relative abundance for deuterium, critical 
in fixing the baryon fraction, of D/H = (3.0 ± 0.4) • 
lO^'^ jlj]. Recently a new measurement of deuterium 
in the damped Lyman-a system was presented p5| , 
leading to the weighted abundance D/H = (2.2 ± 
0.2) • 10^^. We use the value from [Q in our analysis; 
the use of leads to even stronger bound. 

In the standard BBN scenario, the primordial abun- 
dances are a function of the baryon density rj ^ Q,bh? 
only. To constrain the energy density of a primor- 
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FIG. 2. 1,2 and 3cr likelihood contours in the 
(Jlfc/i^, r2^(lMeV)) parameter space derived from '^He and 
D abundances. 



dial field a T - MeV, we modified the standard BBN 
code including the quintessence energy compo- 
nent rj^. We then performed a likelihood analysis in 
the parameter space [VL^h? , fi^^^) using the observed 
abundances Yjjg and D/H. In Fig. ^we plot the 1,2 
and 3tT likelihood contours in the (ri{,/i^, Q.^^^) plane. 

Our main result is that the experimental data for 
^He and D does not favour the presence of a dark 
energy component, providing the strong constraint 
ri0(MeV) < 0.045 at 2cf (corresponding to A > 9 for 
the exponential potential scenario], strengthening sig- 
nificantly the previous limit of ||l^, ^^(JsAeV) < 0.2. 
The reason for the difference is due to the improve- 
ment in the measurements of the observed abun- 
dances, especially for the deuterium, which now corre- 
sponds to approximately ATVpff < 0.2 — 0.3 additional 
effective neutrinos (see, e.g. [18|), whereas Ref. Q 
used the conservative value ATV^g- < 1.5. 

One could worry about the effect of any underesti- 
mated systematic errors, and we therefore multiplied 
the error-bars of the observed abundances by a factor 
of 2. Even taking this into account, there is still a 
strong constraint il^(MeV) < 0.09 (A > 6.5) at 2a. 

Constraints from CMB. The effects of a scaling 
field on the angular power spectrum of the CMB 
anisotropics are several . Firstly, if the energy den- 
sity in the scaling quintessence is significant during 
the radiation epoch, this would change the equality 
redshift and modify the structure of the peaks in the 
CMB spectrum (see e.g. [|9|). 

Secondly, since the inclusion of a scaling field 
changes the overall content in matter and energy, the 
angular diameter distance of the acoustic horizon size 
at recombination will change. This would result in a 
shift of the peak positions on the angular spectrum. 
It is important to note that this effect does not quali- 
tatively add any new features additional to those pro- 
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FIG. 3. CMB anisotropy power spectra for the Al- 
brecht-Skordis models with A = 3, (jAo = 87.22 and 
A = 0.009 (long dash), A = 8, <^o = 32.329 and A = 0.01 
(short dash) (both with Q.^ = 0.65 and h — 0.65), and 
a cosmological constant with ~ 0.65, N,^ — 3.04 (full 
line) and iV,^=7.8(dash-dot). 
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FIG. 4. Matter power spectra for 3 models in fig. 3. 
The predictions support those in the CMB spectra, the 
quintessence model in agreement with BEN A = 8 (short 
dash) mimics the Ha = 0.65 spectrum with A^^=3.0 (full 
line). The model with A = 3 (long dash) is in clear dis- 
agreement with observations. 



duced by the presence of a cosmological constant Ig^l • 
Third, the time-varying Newtonian potential after 
decoupling will produce anisotropics at large angu- 
lar scales through the Integrated Sachs- Wolfe (ISW) 
effect. Again, this effect will be difficult to disentan- 
gle from the same effect generated by a cosmologi- 
cal constant, especially in view of the affect of cos- 
mic variance and/or gravity waves on the large scale 
anisotropics. 

Finally, the perturbations in the scaling field about 
the homogeneous solution will also slightly affect the 
baryon-photon fluid modifying the structure of the 
spectral peaks. However, this effect is generally neg- 
ligible. 

From these considerations, supported also by re- 
cent CMB analysis pl| , p2[ , we can conclude that the 
CMB anisotropics alone cannot give significant con- 
straints on Wtfi at late times. If, however, is signif- 
icant during the radiation dominated epoch it would 
leave a characteristic imprint on the CMB spectrum. 
The CMB anisotropics can then provide a useful cross 
check on the bounds obtained from BBN. 

To obtain an upper bound on fi^ at last scatter- 
ing, we perform a likelihood analysis on the recent 
Boomerang and DASI |^ data. The anisotropy 
power spectrum from BOOMERang and DASI was 
estimated in 19 bins between £ ^ 75 and £ — 1025 
and in 9 bins, from £ = 100 to ^ = 864 respec- 
tively. Our database of models is sampled as in |^ , 
we include the effect of the beam uncertainties for the 
BOOMERanG data and we use the public available 
covariance matrix and window functions for the DASI 
experiment. There are naturally degeneracies between 
fim and JIa which are broken by the inclusion of SNla 
data 10]. 



By finding the remaining "nuisance" parameters 
which maximise the likelihood, we obtain Qcf, < 0.39 
at 2(7 level during the radiation dominated epoch. 
Therefore, while there is no evidence from the CMB 
anisotropics for a presence of a scaling field in the ra- 
diation dominated regime, the bounds obtained are 
actually larger than those from BBN. 

In Fig. 3 we plot the CMB power spectra for 4 alter- 
native scenarios. The CMB spectrum for the model 
which satisfies the BBN constraints is practically in- 
distinguishable from the spectrum obtained with a 
cosmological constant. Nonetheless, if the dark energy 
component during radiation is significant, the change 
in the redshift of equality leaves a characteristic im- 
print in the CMB spectrum, breaking the geometri- 
cal degeneracy. This is also found when considering 
non- minimally coupled scalar fields [Q, even when 
the scalar is a small fraction of the energy density at 
last scattering. In the minimally coupled models con- 
sidered here, this is equivalent to an increase in the 
neutrino effective number, i.e. altering the number of 
relativistic degrees of freedom at last scattering. 

In Fig. 4 we have plotted the corresponding mat- 
ter power spectra together with the decorrelated data 
points of Ref. ||2^. As one can see, the model with 
A = 3 is in disagreement with the data, producing 
less power than the model with A = 8, with this last 
one still mimicking a cosmological constant. The less 
power can be still explained by the increment in the 
radiation energy component which shift the equality 
at late time and the position of the turn-around in the 
matter spectrum towards larger scales. A bias fac- 
tor could in principle solve the discrepancy between 
the A = 3 model and the data, however, the mat- 
ter fluctuations over a sphere of size Sh^^Mpc are 
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(Tg ~ 0.5 to be compared with the observed value 
cTg 0.56f7„''-47 ^ 0.9 [!§. Weak lensing observa- 
tions | |27[ | may open up further opportunities to con- 
strain quintessence models even more tightly through 
the matter power spectrum. 

Conclusions. We have examined BBN abundances 
and CMB anisotropics in a cosmological scenario with 
a scaling field. We have quantitatively discussed 
how large values of the fractional density in the scal- 
ing field 17^ at T ~ IMeV can be in agreement 
with the observed values of "^He and D, assuming 
standard Big Bang Nucleosynthesis. The 2a limit 
V,^{lMeV) < 0.05 severely constrains a wide class of 
quintessential scenarios, like those based on an expo- 
nential potential. For example, for the pure expo- 
nential potential the total energy today is restricted 
to = |r20(lMeV) < 0.04. Our 2a limit on the 
A parameter, could also place useful constraints on 
other dark energy models. In the case of the Albrecht- 
Skordis model, for example, combining our result with 
the condition AX^ < 1, one finds that A < 0.02 in or- 
der to have an eternal acceleration. Furthermore, if we 
want to have Sl^ > 0.65, then one must have 0o < 29. 

As mentioned earlier, our BBN constraint is limited 
to models assuming standard Big Bang Nucleosynthe- 
sis, where the scaling field simply adds energy density 
to the expanding universe. 

The bound on ^^(MeV) can be also weakened by 
introducing new physics which might change the elec- 
tron neutrino distribution function. Such distortions 
will alter the neutron-proton reactions, and subse- 
quently the final abundances. Such new physics could 
take the form of heavy decaying neutrinos, or light 
electron neutrinos oscillating with a sterile species, 
both of which can lead to fewer effective degrees of 
freedom (see e.g. [ p8|j29| and references therein). 

There are several quintessence models which evade 
the BBN bound, and let us mention a few. The sim- 
plest way is to modify the standard model of reheating 
in order to have late-entry of the field in the attractor 
solution, after BBN In the 'tracking oscillating 

energy' of Ref. |Q one can chose parameters which 
both let fi^ be small at BBN times and big today. 
However, since the probability of the randomly se- 
lected parameters decreases rapidly for stronger BBN 
constraints, this model may not appear too natural (in 
the language of ) . Another class of models which 
evade the BBN bound are the models exhibiting kina- 
tion at early times. This means that Qct> is suppressed 
at early times, taking a value well below that made 
by the BBN constraint . 

All these models are compatible with our la con- 
straint obtained from CMB data but they nonethe- 
less leave a characteristic imprint on CMB and also 
large scale structure. Inevitably therefore future data 
from satellite experiments like MAP or Planck, and 
measurements of the matter power spectrum using 
weak lensing and the Digital SLOAN survey will en- 
able even tighter limits to be placed on the presence 



of a scalar quintessence field in the early universe. 
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